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The phase diagram of the Kondo lattice Hamiltonian with ferromagnetic Hund's coupling in the 
limit where the spin of the localized t2g electrons is classical is analyzed in one dimension as a 
function of temperature, electronic density, and a direct antiferromagnetic coupling J' between the 
localized spins. Studying static and dynamical properties, a behavior that qualitatively resembles 
experimental results for manganites occurs for J' smaller than 0.11 in units of the Cg hopping 
amplitude. In particular a coexistence of ferromagnetic and antiferromagnetic excitations is observed 
at low-hole density in agreement with neutron scattering experiments on La2-2xSri+2xMn207 with 
X — 0.4. This effect is caused by the recently reported tendency to phase separation between hole- 
rich ferromagnetic and hole-undoped antiferromagnetic domains in electronic models for manganites. 
As J' increases metal-insulator transitions are detected by monitoring the optical conductivity and 
the density of states. The magnetic correlations reveal the existence of spiral phases without long- 
range order but with fairly large correlation lengths. Indications of charge ordering effects appear 
in the analysis of charge correlations. 
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I. INTRODUCTION 

The study of manganese oxides with perovskite-type 
structure such as Lai-xCaxMnOa and Lai-xSrxMnOa 
has received new emphasis since the discovery of colossal 
magnetoresistance in these compounds ^ . Experiments 
have revealed a rich phase diagram with antiferromag- 
netic insulating, ferromagnetic metallic, and charge or- 
dered regions P,^. Recently, layered compounds with 
properties similar to the 3D materials have been syn- 
thesized: the new compounds are described by the gen- 
eral formula (A, B)n+iMnn03n+i, where n is the num- 
ber of Mn — O planes per unit cell (n — oo corresponds 
to the 3D case) The appearance of ferromagnetism 
upon hole doping away from a: = has been attributed 
to the double-exchange (DE) mechanism between the Cg 
and t2g electrons [gj-Q. However, this model is not suf- 
ficient to describe the rich phase diagram observed ex- 
perimentally. The large regions in the phase diagram 
which present charge ordering are not contained in the 
DE model. Features such as the antiferromagnetic phase 
in the low eg-electronic density limit (a; ^ 1), or the re- 
lated experimentally observed antiferromagnetic ground 
state detected in Lai-xCaxMnOa for x > 0.55 are diffi- 
cult to explain by DE ideas. The Neel temperature of this 
state reaches a maximum value of 260K and decreases to 
120K at X = 1 |]. 

The main goal of this paper will be to study the phase 
diagram of a modified ferromagnetic (FM) Kondo lat- 
tice Hamiltonian in order to understand qualitatively 
what properties of the manganites are solely due to elec- 
tronic interactions. An analysis including Jahn- Teller 
phonons |0] will be postponed for a future publication. To 



achieve our goal the ferromagnetic Kondo lattice Hamil- 
tonian will be here modified by introducing a direct anti- 
ferromagnetic exchange interaction, J', between the spins 
of the localized t2g electrons ||9|,|l0||. This coupling will 
certainly produce an antiferromagnetic phase at x — I, 
as experimentally observed, which does not appear in 
the standard FM Kondo model. The addition of this 
coupling is certainly not unphysical since it may origi- 
nate from small hopping amplitudes for the t2g electrons. 
Experimentally such an effective antiferromagnetic ex- 
change interaction J' has been studied in the 2D com- 
pounds and 2J'S « lOOK has been estimated 

Among the main issues that will be addressed in the 
present manuscript is the influence of the direct cou- 
pling J' among the t2g electrons on the recently reported 
regime of phase separation detected in the ferromagnetic 
Kondo model This effect, as it will be discussed 

later, could be a key feature to understand some of the 
recent neutron scattering results for the manganites. 

The organization of the paper is the following: in Sec- 
tion II the Hamiltonian will be introduced. Phase sepa- 
ration will be studied in Section III. Magnetic and charge 
ordering will be discussed in Section IV. In Section V the 
existence of a Fermi surface is studied by monitoring n{k) 
and the spectral function A{k, u). The optical conductiv- 
ity, the Drude weight and metal-insulator transitions are 
presented in Section VI. In Section VII the phase diagram 
is proposed. Section VIII is devoted to summarizing the 
results. 
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II. THE MODEL 



III. PHASE SEPARATION 



The parent compound LaMnOs contains Mn'^+ ions in 
a iigSg configuration. The electrons may be viewed as 
a locahzed S* = 3/2 spin, while the electrons are mo- 
bile. Here one single orbital for the Cg electrons will be 
considered. To describe the doped manganites a Kondo 
lattice model with ferromagnetic Hund coupling has been 
proposed In spite of its apparent simplicity, this 

model is still difficult to study accurately. In order to 
simplify its analysis the localized S — 3/2 spins will be 
replaced by classical spins with | | 1 Q . In addition, 
an antiferromagnetic coupling J' will be added between 
the classical spins to obtain the antiferromagnetic phase 
at a; = 1 as explained in the Introduction PJlC|]. Fi- 
nally, since the standard FM Kondo Hamiltonian in 3D 
is difficult to study numerically due to its large number 
of degrees of freedom, some other simplifications must 
be introduced. Since previous studies have convincingly 
shown that the basic features of the phase diagram in 
three dimensional problems appear in two and one di- 
mensions as well llQ], results in one dimension will be 
here presented. A careful study of the ID problem in 
chains with up to 40 sites will be carried out using Monte 
Carlo (MC) techniques. 

Then, the Hamiltonian studied in this paper is given 

by 



-Jh c\ a(Ja,pCi^p ■ S 



h.C.) 



(1) 



where c| ^ creates an Cg electron at site i with spin pro- 
jection a, a are Pauli matrices, Si is the total spin of the 
t2g electrons assumed localized and classical, the sum (ij) 
runs over pairs of nearest-neighbor lattice sites, t is the 
nearest-neighbor hopping amplitude for the Cg electrons, 
Jh > is the Hund coupling, and J' > is an anti- 
ferromagnetic coupling between the localized spins. The 
density (n) — I ~ x of itinerant Cg-electrons is controlled 
by a chemical potential fi. In the following t — 1 will 
be used as the unit of energy and Jh will be fixed at 
8 since phenomenologically Jh ^ t. Periodic boundary 
conditions (PBC) will be used. The model will be stud- 
ied with the Monte Carlo technique already described in 
Ref. This method basically amounts to integrating 
exactly the eg electrons using library subroutines (since 
it corresponds to a one-electron problem in the back- 
ground of the localized spins), supplemented by a Monte 
Carlo simulation on the t2g classical spins. The tech- 
nique does not have "sign" problems and it allows to 
reach low temperatures at any density. The neglect of 
the direct Coulomb interactions among the Cg electrons 
can be justified recalling that a large Hund coupling ef- 
fectively reduces double occupancy, similarly as a large 
on-site Hubbard repulsion does. 



In Ref. ||T^ the phase diagram of model Eq.(l) has been 
studied for the special case J' = 0. In that paper it was 
found that at low temperatures the antiferromagnetic in- 
sulating phase that appears at a; = is phase separated 
from a metallic ferromagnetic phase that becomes stable 
for X > 0.25. Here the effect that a finite J' has on this 
important result will be investigated. 

To study phase separation the density (n) as a function 
of the chemical potential fi has to be monitored. Working 
with a chain of L sites and at zero temperature (T), the 
density (n) can take only the values ^ with m running 
from to 2L. However, at finite T the effects of thermal 
fluctuations allow for intermediate values of (n) to exist, 
and this complicates the study of phase separation since 
sharp discontinuities are usually transformed into rapid 
crossovers. This behavior is clearly observed in long MC 
runs when tunneling between several of the allowed states 
occurs. In particular during MC runs of 40,000 iterations 
at low temperature the system was observed to tunnel 
several times between states with a sharply defined value 
of (n) for particular values of the chemical potential. In 
Fig.l results for a Monte Carlo run for /x = —6.97 on an 
L — 20 chain and at T = 1/75 are presented as an il- 
lustration. The density clearly fiuctuates between states 
with 15, 16, 17, 18 and 19 electrons. This particular 
behavior allowed us to measure ground state properties 
for most of the possible filling levels by simply carrying 
out measurements of, e.g., the density and energy in the 
plateaus observed in Fig.l. With this information (i.e. 
energy vs density without temperature contaminations) 
it was possible to perform an accurate analytical calcula- 
tion of (n) vs /i at T = using the Maxwell construction. 
This method allows us to obtain ground state informa- 
tion from a finite temperature simulation. 
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FIG. 1. The total number of itinerant electrons Ne as a 
function of the Monte Carlo steps for an L = 20 chain with 
Jh = 8, J' = 0.05, T = 1/75 and = -6.97. 
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In Fig. 2a the relation (n) vs /i at T = obtained ana- 
lytically from the MC results at J' = 0.05 on an L = 20 
lattice (solid line) is presented, and it is compared with 
the results obtained directly from Monte Carlo simula- 
tions at T = 0.01 (open circles) The T = curve 
shows that one regime of phase separation (denoted by 
PSl) occurs between regions with (n) — 1 (i.e., x = 0) 
and (n) = 0.65. This behavior is similar to the results re- 
cently reported in Ref. |Q working at J' = 0. However, 
the same figure shows the novel result that there is a sec- 
ond region of phase separation (denoted by PS2) in the 
limit of low e^-electronic density between the empty sys- 
tem (n) = 0.0 and an electron-rich phase with (n) = 0.25. 
Through a careful analysis it has been found that phase 
separation at low density also occurs for smaller values 
of J' such as 0.02. 




FIG. 2. The density of itinerant electrons (n) as a function 
of the chemical potential ^ on L = 20 chains with Jh ~ 8 
and (a) J' = 0.05, (b) J' = 0.02, (c) J' = 0.1, (d) J' = 0.15 
and (e) J' = 0.25. The open circles indicate Monte Carlo 
results at T = 0.01. The lines are obtained using the Maxwell 
construction and they indicate T = behavior. 



It is important to remark that at J' = indications 
of a rapid crossover at low densities exist but this effect 
does not reach the full character of phase separation since 
at X = 1 the t2g spins are non-interacting in this special 
case. In Fig.2b-e, (n) vs /i curves at T = for different 
values of J' using chains with L = 20 sites are presented 
(obtained using the methods discussed above that al- 
lowed us to remove the temperature contaminations from 
the results). Additional MC runs using open boundary 
conditions and different lattice sizes close to the phase 
separation critical chemical potentials have shown that 
the results are not affected by strong finite size effects, 
and thus we are confident that they already represent 
the bulk limit. In general, it was observed (Fig. 2) that 



as J' increases the range of stable densities decreases, 
and for J' > 0.25 the only stable density besides (n) = 1 
and is (n) = 0.5. Evenmore, for values of J' larger 
than 0.6 only the empty and the half-filled phases are 
stable. This remarkably strong tendency to phase sepa- 
ration can be qualitatively understood in the limit where 
J' is larger than t, both still smaller than the Hund cou- 
pling. In this regime the antiferromagnetic order caused 
by J' among the localized spins induces a similar an- 
tiferromagnetic tendency among the itinerant electrons 
(again, note that the Hund coupling is assumed to be 
the largest scale). Then, in the conduction band an ef- 
fective one band t — J' model with J' /t ^ 1 dominates 
the physics. It is known that in this regime the t — J' 
model phase separates between an empty and half-filled 
phases in all dimensions | [l6[ |. This can also be under- 
stood easily by considering the atomic limit of Eq.(l), 
i.e., when t Jh , J' ■ In this case the lowest energy of 
a state with Ne Cg electrons, including the chemical po- 
tential, is given by —JmN^ + J'L — fiN^ ii < L and 
- Jh(2L - Ne) + J'L - ^j,Ne if Ne > L. The system is 
half- filled for /i = 0. When /i = — Jh (Jh) the energies of 
all the levels with N^ < L {N^ > L) become degenerate 
vfiihE = J'L {{J'-2Jh)L) . For^ < - Jh {p > Jh) the 
state with the lowest energy E — J'L {{J' — 2Jh)L) is the 
one with N^ = {2L). For — Jh < fJ, < Jh the ground 
state is the state with N^ — L. Thus, at /i = — Jh (Jh) 
phase separation between the half-filled and the empty 
(fully doubly occupied) state occurs. Similar tendencies 
have been recently observed in one band models for man- 
ganites Q. 

In Fig. 3 the stable densities as a function of J' are 
presented. A very similar result was obtained using open 
boundary conditions. In the following sections the prop- 
erties of the stable regions in the phase diagram will be 
studied. 
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FIG. 3. The stable and unstable regions in the J'-{n) 
plane are indicated for Jh — 8. 
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IV. MAGNETIC AND CHARGE ORDERING 



A. Magnetic 

Some of the most important properties of the man- 
ganite materials are related to the magnetic order of the 
electrons since a remarkable interplay between ferromag- 
netic and antiferromagnetic phases has been observed in 
the experiments Actually recent neutron-scattering 
measurements performed on La2-2xSri_|_2xMn207 with 
X = 0.4 have reported the coexistence of short-range 
antiferromagnetic correlations with strong ferromagnetic 
fluctuations p| . In particular, antiferromagnetic fluctu- 
ations one order of magnitude smaller in intensity than 
the ferromagnetic ones were observed between lOOK and 
200K, i.e., above the Curie temperature. The antifer- 
romagnetic fluctuations disappear when ferromagnetic 
long-range order is established. 

Motivated by these experiments, in the present section 
we will study the behavior of the spin-spin correlation 
functions among the classical spins defined as 



1 



(2) 



where the notation is standard. The structure factor S{q) 
is given by the Fourier transform of the spin-spin corre- 
lations: 



'^c.(r). 



(3) 



where the momentum q can take the values iim/L with 
n nmning from to i — 1. In Fig. 4, S{q) aX T — 1/75, 
J' = 0.05 and for densities (n) — 0, 0.35, 0.5, 0.65 and 
1 is presented. It is clearly observed in the figure that 
for (n) = and 1 the structure factor peaks at momen- 
tum q = n denoting antiferromagnetism in the system 
in agreement with the experimental behavior This 
is certainly the effect that we wanted to induce by the 
inclusion of a J' coupling, namely the existence of AF 
correlations at both ends of the eg density range rather 
than only at x = 0. At the remaining stable values of 
(n) the peak occurs a.t q — or q = 27r/20 signaling the 
presence of ferromagnetism. We observed that when the 
number of itinerant electrons is even and periodic bound- 
ary conditions are used, a kink appears in the ground 
state separating two ferromagnetic regions with opposite 
spin, jlql This kink causes the peak in S{q) to move to 
2tt/L [|19[| . As a result of this effect a kink appears at 
quarter filling when the length of the chain L is a mul- 
tiple of 4. When this happens we have verified that the 
peak is at g = for other values of L and also using 
open and antiperiodic boundary conditions. Then, there 
is no doubt that ferromagnetic correlations are dominant 
when the peak in S{q) occurs at or close to q ~ 0. 
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FIG. 4. Structure factor, S{q), as a function of q/n for an 
L = 20 chain at J' = 0.05 and T = 1/75. The densities are 
indicated in the figure. 
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FIG. 5. (a) Structure factor as a function of qj-n on an 
L = 20 chain at T = 1/75, = -6.99 and J' = 0.05. Maxima 
at g = and tt are clearly observed; (b) same as (a) at dif- 
ferent temperatures indicated in the figure; (c) the structure 
factor at g = as a function of the inverse temperature /3 for 
ii — —6.99; (d) same as (c) for g = tt. 



Note that for all the densities studied in Fig. 4 the 
structure factor has a maximum at a single value of 
the momentum indicating that only one kind of mag- 
netic fiuctuations prevail at a given density. However, 
finite temperature measurements at other densities such 
as (n) = 0.75 which corresponds to an unstable density 
in the limit of zero temperature indicate otherwise. As 
it can be observed in Fig. 5a, strong ferromagnetic corre- 
lations coexist with weak antiferromagnetic ones at this 
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density and T = 1/75. This result has similarities with 
the experimental data jl^ discussed before, and we be- 
lieve it is due to the coexistence of hole-undoped antifer- 
romagnetic and hole-rich ferromagnetic domains rather 
than to a ferromagnetic spin polaron in an antiferro- 
magnetic background, as suggested in Ref. |^^. To ex- 
plore this scenario the chemical potential was fixed at 
/X = —6.99 and S{q) was measured at T = 1/25, 1/50, 
1/75, 1/100 and 1/150. The results are presented in 
Fig. 5b. As the temperature is lowered the intensity of 
the peak at g = increases (Fig. 5c) while the peak at 
q = n decreases (Fig.5d). The density reaches the value 
0.65 at the lowest temperature indicating that the sys- 
tem at 0.75 is unstable. The qualitative similarities of our 
results with the experimental data of Ref. ||ll[] is clear. 
The relative intensity of the ferromagnetic and antiferro- 
magnetic peaks can be tuned by properly selecting the 
density at finite temperature in the region that will be- 
come unstable as T ^ 0. 
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FIG. 6. Structure factor as a function of the momentum 
for different densities at J' = 0.15 and T = 0.01. 



Ferromagnetism in the low-temperature stable inter- 
mediate densities between the two regimes of phase sepa- 
ration PSl and PS2 is observed up to J' « 0.1. However, 
for larger values of J' weak incommensurate correlations 
were observed. In Fig. 6, S^q) at J' = 0.15 and for densi- 
ties (n) = 0.35, 0.5 and 0.6 is presented. A low-intensity 
peak at (n) = 0.35 and 0.6 suggests a tendency to short- 
range incommensurate order with momentum q ~ 0.37r 
and 0.57r, respectively. Note, however, that at (n) — 0.5 
a strong peak is observed at g = 7r/2. This indicates the 
existence of a spin twisted ground state (a spiral state 
with pitch q = 7r/2). The structure factor at (n) = 0.5 
in chains of sizes ranging from L = 10 to 30 was stud- 
ied and it is presented in Fig. 7a. A large increase of 
the peak intensity with the lattice size was not observed. 
However, in all the chains studied we found that the cor- 
relation length for the twisted order although finite seems 



fairly large as it can be observed in Fig. 7b. The inten- 
sity of the incommensurate correlations at quarter-filling 
reaches a maximum at J' k, 0.15. As J' increases further, 
incommensurability at g = 7r/2 weakens and at J' « 0.35 
the ground state becomes antiferromagnetic, which is the 
expected behavior at the stable densities in the large J' 
limit. 
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FIG. 7. (a) Structure factor as a function of g/vr at 
T = 1/75, J' = 0.15 and (n> = 0.5 for L=10, 16, 18, 20, 
24, 26, 28 and 30. A peak at g = 7r/2 appears; (b) real part 
of the spin correlation function a;(r) multiplied by as a 
function of the distance r for L = 24 (triangles), 26 (circles) 
L = 28 (squares) and 30 (diamonds). 



B. Charge 

The existence of commensurate and incommensurate 
charge ordered phases has been reported in the experi- 
mental literature for Lai-xCaxMnOa ||]. It was observed 
that for values of 0.5 < a; < 0.67 antiferromagnetism co- 
exists with commensurate charge order, and charge in- 
commensurability was detected in a ferromagnetic phase 
at higher temperature. This behavior is very puzzling 
since ferromagnetism and charge ordering were not ex- 
pected to appear together, particularly based on analysis 
using the double exchange model. In order to analyze the 
effect of the interplay of double exchange and antiferro- 
magnetic interactions in charge ordering for the itinerant 
fcrmions, the Fourier transform N(q) of the charge cor- 
relations was studied. The charge correlation functions 
are given by 



{n)){n 



i-\-r 



(4) 



where rii is the number operator at site i for the itiner- 
ant fermions. Computationally we found that in the an- 
tiferromagnetic phase at half-filling N{q) peaks at g = tt 
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rather than g = 0. This behavior is due to charge corre- 
lations that are negative (charge repulsion) at very short 
distances rather than to long range ordering, and it is 
in agreement with previous calculations In the fer- 
romagnetic phase we observed that N{q) behaves as for 
non-interacting spinless fermions, i.e. it reaches a max- 
imum at q — 2k^J — 2'K{n) [k^J = 7r(n) is the Fermi 
momentum for free spinless fermions) as q grows from 0, 
and then it becomes flat (Fig.8-a). This indicates that 
the value of the Hund coupling, Jh = 8, is sufficiently 
strong for the itinerant spins to be parallel to the local- 
ized ones, behaving effectively as non-interacting spinless 
fermions in the charge channel. 
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FIG. 8. (a) Fourier transform of the charge correlation 
functions as a function of q/n at the densities indicated in 
the figure for J' = 0.05 on a L = 20 chain at T = 1/75; 
(b) Charge correlation function as a function of distance for 
(n) — 1; (c) same as (b) for (n) = 0.65. The small full 
circles are results obtained using Eq.(5); (d) same as (c) for 
(n) — 0.55; (e) same as (c) for (n) = 0.5. 



In Fig.8b-e the charge correlations as a function of 
distance are presented. At half- filling (Fig. 8b) the only 
non-negligible correlations arc the nearest-neighbor ones 
which are negative. This effect induces a maximum at 
g = TT in N(q). In the ferromagnetic region the charge 
correlations show oscillatory behavior. In particular, at 
quarter filling (Fig.8e) the correlations are zero at even 
distances and negative at odd ones. Although the sys- 
tem is metallic, this produces a maximum of N{q) at 
q = TT similar to the one experimentally observed along 
the a axis of Ca-doped manganites with 50% doping [|j . 
The behavior of the charge correlation for non-interacting 
spinless fermions is given by 



n{r) 



-1 



cos(2fc 



sf 



(5) 



This expression is strictly valid in the limit L ^ oo and 
for r ^ but as it can be observed from Fig. 8c, d and e, 



where the correlations obtained from Eq.(5) are denoted 
with filled circles, they are in very good agreement with 
the numerical results for L ~ 20. Thus, apparently this 
is a system that behaves like a ferromagnet in the spin 
channel, but is non-interacting in the charge channel. 

While the qualitative behavior of N{q) at half-filling 
is independent of J' since the ground state is always an- 
tiferromagnetic, a different behavior is observed at in- 
termediate densities when the ground state is no longer 
ferromagnetic. In Fig. 9a, N{q) is presented for J' = 0.15. 
N[q) now reaches a maximum at g = 2A;|/ and afterwards 
the intensity forms a plateau at a lower value. In Fig.9b-e 
the charge correlations as a function of the distance are 
shown. It is clear that the peaks in N{q) are caused by 
oscillatory behavior of the charge correlations. It appears 
that a tendency to short-range charge ordering exists in 
the model for densities close to quarter-filling. 
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FIG. 9. (a) Fourier transform of the charge correlation 
functions as a function of q/n at the densities indicated in 
the figure for J' = 0.15 on a L = 20 chain at T = 0.01; 
(b) Charge correlation function as a function of distance for 
(n) = 1; (c) same as (b) for {n) = 0.65; (d) same as (c) for 
(n) = 0.55; (e) same as (c) for {n) — 0.5. 



V. FERMI SURFACE AND SPECTRAL 
FUNCTIONS 

In previous sections, regions in parameter space with 
a variety of magnetic and charge structures were found. 
Here the dynamical properties of these phases will be 
studied. A calculation of the spectral function A{k,u!) 
and the study of the momentum distribution function 
n(fc) of the Cg electrons will allow us to determine, at 
least qualitatively, whether a Fermi surface exists and, 
thus, whether the phase is insulating or metallic (this 
analysis will continue in the following section where the 
Drude weight of the optical conductivity will be dis- 
cussed). In principle, the spectral function can be coni- 
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pared with angle-resolved photoemission (ARPES) mea- 
surements performed in Lai 2Sri gMn207 and in 
Lai„,SrxMn03 ||. 

The spectral function A{k,uj) is defined as 



1 



ImG{k, w), 



(6) 



where G(fc, lo) is the one-particle Green's function for the 
eg electrons. Since we are performing a Monte Carlo cal- 
culation on the classical spins only, the time-dependent 
Green's function can be straightforwardly calculated in 
real-time since the fermions do not interact among them- 
selves directly, but only with the classical spins. The 
spectral functions obtained with this procedure only con- 
tain small controlled statistical errors. This should be 
contrasted against other calculations that use the usu- 
ally uncontrolled Maximum Entropy technique. 

Studying the density of states for the Hamiltonian de- 
fined in Eq.(l) we observed low-energy and high-energy 
bands. The separation between these two bands is regu- 
lated by the Hund coupling Jh, i.e., the lower (upper)- 
band includes states where the spin of one of the itinerant 
electrons at an arbitrary site i is parallel (antiparallcl) to 
the classical spin at the same site of the lattice. The gap 
is about 2Jh- Each of the two bands has an internal 
structure and a bandwidth which are regulated in part 
by the antiferromagnetic coupling J' . 



A. Half-filling 



In Fig. 10a the density of states N{lo) corresponding to 
J' = 0, half-filling, and with T = 1/100 on a L=20 chain 
is presented, The bandwidth of the upper and lower 
bands is ~ \M and the gap separating them is ^ 14t. 
The chemical potential lies in the gap indicating that 
the system is an insulator. In Fig. 10b, A{k,uj) is shown. 
Both below and above the chemical potential, spectral 
weight is observed for all of the momenta k, rather than 
only for k below or above the Fermi momentum as in 
non-interacting electronic systems. This is reminiscent 
of the antiferromagnetically induced features observed at 
half-filling in studies of the one band Hubbard model for 
cuprates |2^. Below the chemical potential there is a 
sharp peak for k smaller than 7r/2 at the bottom of the 
lower band. The dispersion of this feature is about t. The 
corresponding opposite behavior is observed in the upper 
band. If J' is increased, the bandwidth of the upper and 
lower bands is slightly reduced due to the increase in the 
AF correlations, as it can be observed in Fig. 11 where 
the half- filled density of states is presented for J' = 0, 
0.05, 0.1, 0.15 and 0.25. 



N(co) 




A(k,co: 




FIG. 10. (a) Density of sates N(uj) as a function oi u — ^ 
for Jif = 8, J' = at (n) = 1 for L = 20 at T = 0.01; (b) the 
spectral function A{k, lu) as a function oi u — ^ for the values 
of k available on an L = 20 lattice. The parameters are as in 
(a). 
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FIG. 11. Density of states as a function of lj — /i at 
half-filling on L = 20 chains at T = 0.01 for the values of 
J' indicated in the figure. 



B. Quarter- filling 

The spectral functions change drastically in the ferro- 
magnetic phase. First let us discuss results for density 
(n) = 0.5 (quarter-filled). Jh still regulates the gap be- 
tween the upper and lower bands. We will concentrate 
our attention on the lower band, since symmetrical fea- 
tures are observed in the upper band. In Fig. 12a, A{k, lo) 
aX J' — Q is shown at quarter-filling and T = 1/75 on a 
i = 40 chain (similar results were obtained for shorter 
chains). The chemical potential lies now in the middle of 
the lower band indicating that the system is a metal with 
a Fermi surface at kp — 7r/2. Notice that kp results to be 
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equal to kp . For each momentum k there is a sharp peak 
which is naturally associated with a quasiparticle peak. 
The dispersion resembles results for a system of free spin- 
less fermions (indicated by a dashed line in the figure), 
in agreement with the behavior of the charge correlations 
described in the previous section. The bandwidth is At. 
The density of states can be seen in Fig. 12b. It resem- 
bles N{uj) for spinless fermions at half-filling which for a 
finite length chain typically presents several spikes that 
become a smooth function only for large enough chains. 
Thus, the spikes observed in Fig. 12b near the chemical 
potential are likely finite size effects, and the quasiparti- 
cle peak crosses the chemical potential as in a standard 
metal. Similar behavior is observed for J' up to approx- 
imately 0.11 in the ferromagnetic phase. 



for J' « 0.6 when the quarter-filled state is no longer 
stable. In Fig. 14 the density of states is presented for 
different values of J' at (rt) = 0.5 using L = 20 chains. 
Once again, the spiky structure for J' — 0.1 or smaller 
is a finite size effect that disappears as the size of the 
lattice grows. However, the gap observed for J' > 0.11 is 
not spurious and signals an insulating regime. The open- 
ing of this gap is not just a one dimensional effect. In 
Fig. 15 we present the density of states at quarter-filling 
for Jh = 8 and J' = 0.25 at T = 1/50 on a 4 x 4 and 
a 6 X 6 lattice with PBC. It can be seen that a gap has 
opened in the lower band and the chemical potential lies 
in it, denoting insulating behavior. The finite size effects 
do not seem strong. 



A(k,co; 



A(k,(o) 




N(co) 




FIG. 12. (a) The spectral function A{k, ui) as a function 
of tj — /X for the values of k available on an L = 40 lattice for 
Jh = 8, J' = 0, T = 1/75 at (n) = 0.5. The dashed line is 
the dispersion for spinless fermions; (b) Density of states as a 
function of cij — /i. The parameters are as in (a). 



However, as J' increases and the incommensurate 
phase is reached the behavior of the spectral function 
at quarter-filling changes dramatically. Results for J' = 
0.15 are presented in Fig. 13 (remember that for this value 
of J' strong incommensurate spin correlations were de- 
tected in the study of S{q)). A gap has now opened and 
the chemical potential lies in it indicating that the sys- 
tem became an insulator. As shown in Fig.l3-b the size 
of the gap appears to be independent of the size of the 
chains. Fig.l3-a shows that considerable spectral weight 
appears above and below the chemical potential for mo- 
menta close to the Fermi momentum k^f = 7r/2. The 
dispersion for the features above and below the chemical 
potential in the lower band is of the order of 2.5t. With 
increasing J' the system remains insulating, but it was 
observed that the size of the gap decreases while anti- 
ferromagnetic correlations increase. The gap closes up 




N(co) 




FIG. 13. (a) The spectral function A{k, ui) as a function 
of u) ~ H for the values of k available on an L = 20 lattice for 
Jh = 8, J' = 0.15, T = 1/75 at (n) = 0.5; (b) Density of 
states as a function of ui — fi. The parameters are as in (a) 
but L = 20 and 32. 




-2.0 0.0 2.0 



FIG. 14. Density of states as a function of — /i for dif- 
ferent values of J' at quarter filling for L = 20 and T — 1/75. 
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FIG. 15. Density of states as function of ui — ^ for J' = 0.25 
at quarter filling for 4x4 and 6x6 lattices at T = 1/50 with 
PBC. 
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FIG. 16. (a) The spectral function A{k, uj) as a function 
of a; — ^ for the values of k available on an L = 20 lattice for 
Jh = 8, J' = 0.1, T = 1/100 at (n) = 0.55; (b) same as (a) 
for (n) — 0.45. The dashed line is the dispersion for spinless 
fermions. 



C. Other Densities 



The spectral functions at densities below and above 0.5 
were also studied. In the ferromagnetic region a behav- 
ior resembling free spinless fermions was observed, as it is 
shown in Fig. 16 where A{k,uj) for J' = 0.1 is presented 
at (a) {n) = 0.55 and (b) (n) = 0.45. As pointed out 
in Sec.IV.B, in the charge channel the itinerant electrons 
behave as free spinless fermions with fc|/ = Tr{n). In the 
insulating phase, on the other hand, the density of states 
changes qualitatively away from quarter filling. When iVg 
electrons {N^ holes) are introduced in the quarter filled 
state, spectral weight is redistributed so that a new peak 
develops inside the gap and, thus, a three peak structure 
is observed. The spectral weight in the original two peaks 
changes from L/2 to L/2 — {L/2 — Nh) while the spec- 
tral weight of the new structure is given by 2Ne {2Nh)- 
p6t However, the chemical potential still lies in the gap 
indicating that the system continues to be an insulator 
as it can be seen in Fig. 17 where the density of states 
for J' = 0.15 at different values of (n) is presented. For 
the cases of (n) = 0.55 and 0.45 the spectral function is 
shown in Fig. 18. It can clearly be seen that, for all values 
of the momentum fc, spectral weight is being transferred 
to the middle of the gap. Notice that while the two ex- 
ternal features disperse as a function of the momentum 
k, the peak at the center does not. This third peak is 
located exactly in the middle of the gap and it indicates 
that as electrons or holes are added to the quarter filled 
insulating system, localized states form in the middle of 
the gap. 
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FIG. 17. Density of states N{lu) as a function of u} — fj, 
for Jh — 8, J' = 0.15 at different densities for L = 20 at 
T = 0.01. 
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(0-)j, co-ia, 

FIG. 18. (a) The spectral function A(k,u}) and the density 
of states N{lu) as a function of a; — on an L = 20 lattice for 
Jh = 8, J' = 0.15 at (n) = 0.55 at T = 0.01; (b) same as (a) 
for (n) = 0.45. 



The occupation number n(k) can be calculated directly 
from the Green's functions or just by integrating A(k, uj) 
on UJ up to the chemical potential. In Fig. 19 n{k) for 
quarter filling and different values of J' is presented. A 
sharp change in n{k) suggests the existence of a Fermi 
surface while a smooth behavior is more compatible with 
the opening of a gap. The change of behavior between 
J' — 0.10 and 0.15 is very clear. The analysis of this 
curves complements the previous study of the spectral 
functions. Analogous results were obtained in both cases. 



with the temperature |2^. In particular the appear- 
ance of a pseudogap with increasing temperature was re- 
marked. In Fig. 20 the density of states at quarter-filling 
for J' = is presented at different temperatures. As the 
temperature increases the peaks due to finite size effects 
are smeared out but the density of states retains its quali- 
tative characteristics. There are no indications of a pseu- 
dogap developing as the temperature increases. Thus, it 
appears that some experimental observations may be due 
to effects that go beyond the electronic interactions con- 
sidered in this work. Probably more than one orbital per 
site as well as phonons are needed to account for this 
nontrivial effect. 
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FIG. 20. The density of states as a function of lj — /i at 
different temperatures for Jh = 8, J' = at quarter filling 
on a 40 site chain. 
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FIG. 19. n{k) as a function of k for different values of J' 
for T — 0.01 at quarter filling on an L = 20 chain. 



Some experimental ARPES studies for manganites 
have found a strong dependence of the density of states 
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FIG. 21. The density of states as a function of lj — ^ at 
different temperatures for Jh — 8, J' = 0.12 at quarter filling 
on a 22 site chain. 
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For completeness, the temperature dependence of the 
density of states was also studied at quarter-filling in the 
insulating phase. In Fig. 21, where results for J' = 0.12 
are presented, it can be seen that an increase in temper- 
ature closes the gap in the density of states transforming 
it into a pseudogap. These results suggest that the real 
manganite materials could be located right at the bor- 
der between the metallic and insulating regimes reported 
here. 



VI. OPTICAL CONDUCTIVITY AND DRUDE 
WEIGHT 

A very important characteristic of the manganite ma- 
terials is the large change in resistivity observed exper- 
imentally with small changes in doping, temperature, 
or magnetic fields. In addition, experimental measure- 
ments of the optical conductivity in Lai-xSr^MnOa with 
< X < 0.3 [g^ indicate the existence of a gap in the op- 
tical spectrum at half-filling which becomes a pseudogap 
in the ferromagnetic region at high temperatures. This 
behavior has similarities with those observed in ARPES 
data. To study the extent to which electronic interac- 
tions can account for the optical measurements, here we 
have calculated the optical conductivity defined as: 



^(1-e-^") 



dt 
2^ 



e'"*(>W>(0)), (7) 



where is the current operator in the x-direction given 

by 



3. 



(8) 



with the rest of the notation standard. At w = there is 
a Drude weight with an intensity D given by: 



D 
2^ 



< -r > 

4L 



1 



27re2 



0+ 



dLOIj(uj), 



(9) 



where T is the kinetic energy operator and e is the elec- 
tronic charge. Hereafter we set e = 1. 

In Fig. 22a the Drude weight at density (n) = 0.5 is pre- 
sented as a function of J' for chains of different lengths, 
and at low temperature. The results suggest that the 
Drude weight of the system decreases as J' increases and 
it actually vanishes for J' k, 0.11, i.e., where the metal- 
insulator transition was observed in previous sections by 
monitoring the spectral function and nik). The effect 
of the temperature on the optical conductivity was also 
studied. In Fig. 22b the Drude weight is plotted as a 
function of temperature for several values of the antifer- 
romagnetic coupling J' working at quarter-filling. For 
values of J' that lead to a ferromagnetic ground state at 
low temperature, the Drude weight increases as the tem- 
perature decreases. On the other hand, a decrease in the 



Drude weight with decreasing temperatures is observed 
when the ground state at low temperature is an insulator 
with strong spin-incommensurate correlations. 
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FIG. 22. (a) The Drude weight for (n) = 0.5 as a function 
of J' for chains of different lengths; (b) Drude weight as a 
function of temperature at quarter filling. The values of J' 
are indicated in the figure. 



In Fig. 23 the optical conductivity for J' = at differ- 
ent temperatures is presented. It is clear that the spectral 
weight on the high-energy side, i.e. across the 2Jh gap, 
decreases as the temperature decreases while the spectral 
weight on the low lying band increases. This occurs be- 
cause at high temperature both the spins up and down 
have similar occupations of their lower bands and the op- 
tical gap transitions have contributions from both kind of 
spins. As the temperature is lowered and ferromagnetism 
develops, the density of states for the spins antiparallel 
to the classical spins goes up in energy while the density 
of states for the parallel spins at low energy increases 
providing low energy states for transitions, and as a con- 
sequence the transitions across the 2Jh gap diminish. 
Similar results have been obtained in infinite dimension. 
p8[ This behavior is in qualitative agreement with ex- 
perimental measurements of the optical conductivity in 
doped LaMnOs 
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FIG. 23. Optical conductivity as a function of uj for 
Jh = 8, J' = and quarter filling for different values of the 
temperature; (a) high energy; (b) low energy. A lorentzian 
with width 0.25 was used. 



erant and localized electrons via a Zeeman term, 
-'iEi(Ea,/3 4a(^^)«,/3Ci,/3 + 35'f/2) . A insulator-metal 
transition induced by the magnetic field h was observed 
for values of J' where the ground state is an insulator, as 
it can be seen in Fig. 24. 



VII. PHASE DIAGRAM 

The information about the Kondo Hamiltonian pre- 
sented in the previous sections allows us to obtain the 
low temperature phase diagram as a function of the anti- 
ferromagnetic coupling J' and the density (n) in the limit 
of large Hund coupling. The full phase diagram is pre- 
sented in Fig. 25. For J' smaller than ^0.11 the groimd 
state is a ferromagnetic metal for (n) between ~ 0.25 and 
~ 0.65 and an antiferromagnetic insulator for the empty 
and half-filled systems that are phase separated from the 
intermediate densities. This is the region of the phase 
diagram which is relevant for comparisons with exper- 
imental results for the manganites. When J' becomes 
larger than ~ 0.11 a metal-insulator transition occurs 
and the system becomes an incommensurate insulator at 
intermediate fillings. In particular the ground state at 
quarter-filling appears to be a spin spiral state with a 
maximum at q = 7r/2 in the magnetic structure factor, 
and finite but very long correlation lengths. The allowed 
intermediate densities diminish as J' increases and for 
J' Ri 0.25 only quarter-filling is stable. At J' « 0.35 
the quarter filled ground state becomes antiferromagnetic 
and for J' larger than ~ 0.6 only the half-filled and the 
empty phases are stable. Both of them correspond to 
antiferromagnetic insulators. 
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FIG. 24. Drude weight as a function of magnetic field for 
J' = 0.14, T = 1/75 at quarter filling on an L = 10 chain. 



Finally, for completeness the influence of a mag- 
netic field on our results was also considered. The 
external magnetic field was coupled to the itin- 
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FIG. 25. The phase diagram for the Kondo lattice Hamil- 
tonian with ferromagnetic Hund's coupling Jh = 8. AF in- 
dicates the antiferromagnetic phases, I the insulating region, 
FM the ferromagnetic phase and PS denotes phase separation. 



12 



VIII. SUMMARY 



IX. ACKNOWLEDGEMENTS 



In this paper we have studied a ferromagnetic Kondo 
lattice model with classical localized spins coupled anti- 
ferromagnetically, using Monte Carlo methods. Work- 
ing at large Hund coupling we found a very rich 
phase diagram as the density of itinerant electrons and 
the strength of the antiferromagnetic coupling J' was 
changed. In particular, phase separation between hole- 
undoped antiferromagnetic and hole-rich ferromagnetic 
states was observed at high electronic densities. Recip- 
rocally, phase separation between electron-undoped an- 
tiferromagnetic and electron-rich ferromagnetic regions 
was observed in the limit of low Cg electronic density. At 
intermediate densities, we found a metallic ferromagnetic 
phase for small J' separated from an insulating regime 
with a transition located at J' « 0.11. The insulating 
phase has very strong incommensurate magnetic correla- 
tions with a structure factor that peaks at g = 7r/2 for 
quarter-filling. At J' ~ 0.35 the quarter-filled ground 
state becomes antiferromagnetic and beyond J' « 0.6 
only the half-filled and empty antiferromagnetic phases 
are stable. 

J' = 0.05 is possibly a good value to compare with 
experimental results. We have observed ferromagnetic 
and antiferromagnetic spin correlations coexisting for 
(n) = 0.75 at T = 1/75. As in neutron scattering ex- 
periments the ferromagnetic peak increases as the tem- 
perature decreases, while the antiferromagnetic excita- 
tions vanish in the ferromagnetic phase. In the present 
model this behavior is due to the coexistence of ferro- 
magnetic and antiferromagnetic domains at high tem- 
perature, which are a consequence of the tendency of the 
system to phase separate ]l0| . 

The study of dynamical properties shows the existence 
of interesting metal-insulator transitions as a function of 
temperature and J'. A transference of spectral weight 
from higher to lower energies with decreasing tempera- 
ture, in agreement with experimental results, is observed 
in the optical conductivity. The metal-insulator tran- 
sition with increasing J' is a feature of the model also 
present in two dimensions and it is very likely that the 
transition will also occur in three dimensions. 

Summarizing, the present results show that several 
properties of the manganites, such as the change in 
magnetic ordering as a function of concentration and 
transference of spectral weight from high to low ener- 
gies in the optical conductivity with temperature, may be 
attributed principally to simple electronic interactions. 
However, additional factors will have to be considered to 
explain, for example, the opening of gaps in the density 
of states with increasing temperature. The addition of 
phononic degrees of freedom will probably be relevant 
and we will take them into consideration in future work. 
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